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Elliptic Flow from a Beam Energy Scan: 
a signature of a phase transition to the Quark-Gluon Plasma 
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We employ a relativistic transport theory to describe the fireball expansion of the matter created 
in ultra-relativistic heavy-ion collisions (uRHICs). Developing an approach to fix locally the shear 
viscosity to entropy density 77/s, we study the impact of a temperature dependent rj/s{T) on the 
build-up of the elliptic flow, 112, a measure of the angular anisotropy in the particle production. 
Beam Energy Scan from y^ijvjv = 62.4GeV at RHIC up to 2.76 TeV at LHC has shown that the 
V2{pt) as a function of the transverse momentum pt appears to be nearly invariant with energy. 
We show that such a surprising behavior is determined by a rise and fall of rj/s{T) with a minimum 
at T ~ Tc, as one would expect if the matter undergoes a phase transition or a cross-over. This 
provides an evidence for phase transition occurring in the uRHIC's and a first constraint on the 
temperature dependence of 77/s. In particular, a constant 77/s at all temperatures or a too strong 
T-dependence would cause a breaking of the scaling of V2{pt) with the energy. 
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The main motivation for the ultra-relativistic heavy- 
ion collisions program was to create a transient state of 
a quark-gluon plasma matter (QGP) [l|, 0]. One of the 
main discoveries at RHIC was that such a matter has a 
very low shear viscosity to entropy density ratio rj/s (3| 
close the conjectured lower bound for a strongly interact- 
ing system in the limit of infinite coupling, 77/s — l/47r 
[j|. The first and principal observable indicating such a 
low viscosity is the so called elliptic flow, a measure of the 
anisotropy in the angular distribution of the emitted par- 
ticle defined as V2 = {cos{2ipp)) = {{pi -pi)/ {pi +Py)) 
with ifp being the azimuthal angle in the transverse plane 
and the average meant over the particle distribution. 
When analyzed quantitatively by means of hydrodynam- 
ical simulation [5|, it has been found that the amount 
of V2{pt) observed is nearly consistent with the one of 
a perfect fluid, while an advanced analysis by means of 
viscous hydrodynamics [g, |7| or transport kinetic theory 
[a-[l3| both confirm that the data on V2 at RHIC and 
LHC are consistent with an average Atttj/ s ~ 1 — 3. More 
recently the possibility to measure event- by- event the an- 
gular distribution of emitted particle has made possible 
the measurement of higher harmonics w„ — {cos{rnpp)) 
with n > 2 showing a fast decrease of harmonics with 
n > 3 again cornpatible with a finite but not too large 
value of ?7/s [11, [!l. 

A low value of r]/ s ~ 0.1 in itself is not a direct sig- 
nature of the creation of a QGP and indeed a QGP as 
expected from asymptotic freedom should have an 77/s 
about one order of magnitude larger [16]. It is known 
from atomic and molecular physics that a minimum in 
77/s is expected close to the transition temperature as em- 
phasized in the context of QGP in Refs. [i3jl3- Such a 
minimum can be relatively smooth if one considers a sys- 
tem above the critical point or more pronounced at and 
below the critical point with the possibility to have dis- 



continuity pTHigj . For this reason not an average value 
for 77/s, but rather a phenomenological estimate of its 
temperature dependence, is desired to find a confirmation 
that the matter created undergoes a phase transition. 

From both chiral perturbation theory for a meson gas 
[20I [2l| as well as a transport analysis in uRQMD |22[ 
and the phenomenological analysis performed for heavy- 
ion collisions at intermediate energy (HIC-IE) [23, llj] 
indicate a quite high value Anrj/s > 6 for hadronic matter 
at a temperature T < Tc = 165 MeV, see triangles and 
diamonds in Fig. [TJ At higher T, first data on lattice 
QCD [Hi are compatible with 47r77/s -- 1 - 2 at T -- T^ 
[2a | , see up-triangles in Fig. [U even if uncertainties are 
too large and the calculations have been performed only 
in the quenched approximation, i.e. only in the limit of 
infinite quark masses. 

It has been emphasized by both the STAR Collabora- 
tion at RHIC [li] and the ALICE at LHC [13, [11 (in 
agreement with the measurement done also by CMS (29| 
and ATLAS [30|) that surprisingly the V2{pt) appears 
to be invariant in the very wide colliding energy range 
of 62.4 GeV < ^sii7 < 2.76 TeV. Such an observation 
appears to be quite surprising at first sight because one 
would expect that lowering the energy the contribution 
of hadronic matter over the evolution of the expanding 
matter plays an increasing role in damping the V2{pt), as 
indeed observed for the momentum averaged {V2) [U, 1^- 
It therefore appears a key question to answer the reason 
for the invariance of V2{pt)- 

In this Letter, we show that the invariance of V2{pt) in 
the energy range 62.4 GeV < ^snn < 2.76 TeV is caused 
by a fall and rise of the rj/ s{T) as one would expect if 
the created matter undergoes a phase transition. We em- 
ploy a transport kinetic theory approach, developed to 
perform realistic simulations of HICs keeping the local 
•q/ s{T), to analyze the impact of different temperature 



dependences assumed for the expanding matter. Our 
study reveals that the invariance of W2(pt) at varying 
coUiding energies means that the t]/s{T) has a typical 
" U" shape with a decreasing behavior from the hadronic 
matter and a not too steep rise with temperature in the 
QGP. 

Transport at fixed rj/ s - We have developed in the 
recent years a Relativistic Boltzmann Transport (RBT) 
approach that, instead of focusing on specific microscopic 
calculations or modelings for the scattering matrix, fixes 
the cross section in order to have the wanted rj/ s. This is 
not the usual approach to transport theory that is gener- 
ally employed by starting from cross sections and mean 
fields derived in microscopic models. The motivation for 
our approach is inspired by the success of the hydrody- 
namical approach that has shown the key role played by 
the "q/ s. Therefore on one hand we use the RBT equation 
as an approach converging to hydrodynamics for small 
scattering relaxation time t ^ ap (small rj/s). On the 
other hand the RBT equation is naturally valid also at 
large 77/3 or pr » T (explored in the present work) in 
contrast to hydrodynamics, and avoids uncertainties in 
the determination of the viscous correction, i5/, to the 
distribution function f(x,p), that usually becomes quite 
large at pr > 1.5 GeV j3]| . 

To study the expansion dynamics with a certain 
r]/s{T), we determine locally in space and time the 
total cross section atot according to the Chapmann- 
Enskog theory. For a pQCD inspired cross section. 



da/dt 
approaches 



typically used in parton cascade 
4361 1 ■ this gives: 



V/S = — (P> ^r, = TT-n 

15 15 g{a)atotp 



(1) 



where a = m£)/2T, with mD being the screening mass 
regulating the angular dependence of the cross section 
atot, while g{a) is the proper function accounting for the 
pertinent relaxation time tI'^ — g{a)atotP associated to 
the shear transport coefficient and given by: 
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with Kn-s being the Bessel functions and the function 
h is relating the transport cross section to the total one 
(Ttr{s) = atot h{ml/s) and h{0 = 4C(1+C) [(2C+l)Zn(l+ 

1/0-2]. ^ 

The maximum value of g, namely g{mD —> 00) — 2/3, 
is reached for isotropic cross section and Eq.(IT|) reduces 
to the relaxation time approximation with t^^ — ^~ — 
We have shown in Ref. 



(^trP- 



37| that Eq. 



— 'tr — 

J) correctly 



describes the rj/s of the system in the range of interest 
and it is in good agreement with the Green-Kubo for- 
mula. We notice that in the regime where viscous hydro- 
dynamics applies the specific microscopic details of the 
cross section are irrelevant, and ours is the only effective 
way to employ transport theory to simulate a fluid at a 
given rj/s. 



We solve the RBT equation with the constraint that 
rj/ s{T) is fixed during the dynamics of the collisions in a 
way similar to [3a|, but with an exact local implementa- 
tion as described in detail in '»|. From Eq.(IT]) the cross 
section atot{p,T) determining the wanted value 77/s is 
given by: 
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FIG. 1: Different temperature dependent parametrizations for 
77/s. The orange area takes into account the quasi-particle 
model predictions for rj/s [3^. The different lines indicate 
different T dependencies assumed in the simulation of heavy- 
ion collision. Symbols are as in the legend. See the text for 
more details. 

In our calculation the initial condition is longitudi- 
nal boost invariant flow, but the dynamical evolution is 
3D-I-1. For studying V2 this approximation is adequate, 
although for other collective flow phenomena, like rota- 
tion or turbulence [40, |4l[ more realistic initial conditions 
would be necessary. The initial dN/drj have been chosen 
in order to reproduce the final dNch/dri{b) at mid rapid- 
ity as observed in the experiments at RHIC and LHC 
energies [23, \^- The partons are initially distributed 
according to the Glauber model in coordinate space. In 
the momentum space the distribution is thermal up to 
Pt = 2 GeV and at larger px we include the spectrum 
of non-quenched niinijets according to standard NLO- 
pQCD calculations. In order to fix the maximum tem- 
perature in the center of the fireball, Tmo, we assume 
that it scales with the collision energy according to the 
relation 



1 dN,i, 
tAt drj 



cxT^ 



(4) 



and for the initial time, tq, we ensure that it satisfies the 
uncertainty relation between the initial average thermal 
energy and the initial time by TmoTo ~ 1. Combining 
these two relations one has 



r(Vir) _ dNch/d7]{y^) 
T(^) Y dN,h/dv{^) 



(5) 



as commonly done in hydrodynamical studies [435. Thus 
at 62.4 GeV, 200 GeV and 2.76 TeV the maximum ini- 
tial temperature, Tmo, has the values 290 MeV, 340 MeV 
and 560 MeV respectively. Once the maximum temper- 
ature is fixed, the local temperature profile scales with 
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FIG. 2: Differential elliptic flow V2{pt) at mid rapidity for 
10% — 20% collision centrality. The thin solid line, the dashed 
line and the thick solid line refer to: Au + Au at y/s — 62.4 
GeV and ^/s = 200 GeV and Pb + Pb at ^ = 2.76 TeV, re- 
spectively. On the left panel, the results of the simulation with 
i-KTi/s = 1 during the whole evolution of the system, while on 
the right panel with the inclusion of the kinetic freeze-out are 
shown. See text for more details. The dotted line is the result 
excluding minijets at y^ = 2.76 TeV. 

In the following we will discuss the impact of the tem- 
perature dependence of 77/s, In the cross-over region 
that we identify as the region just below the knee in 
the e/T^ curve [3, i.e. at e < l.SGeV/fm^ and/or 
T < Tq = 1.2Tc, the rj/s should increase linearly at de- 
creasing T matching the estimates from chiral perturba- 
tion theory for a high temperature meson gas [20, l2ll |. 
shown by triangles in Fig. [TJ We notice that the last 
are also comparable to the estimate of 77/s extrapolated 
from heavy-ion collisions at intermediate energies (HIC- 
lE diamonds in Fig. [T]), even if one should consider that 
they refer to a matter with higher baryon chemical po- 
tential fiB- As discussed above, due to the large error 
bars in the IQCD results for ij/s it is not possible to infer 
a clear temperature dependence in the QGP phase. We 
have considered two cases. One with a linear dependence 
Anrj/s = T/Tq (red solid line in Fig. [1]) in a gree ment 
with the indication of IQCD calculation in Ref. [2^ (up- 
triangles) as well as with quasi-particle model prediction 
(orange band) suggesting 77/s ~ T" with a ~ 1 — 1.5 
|39l |45| . The other one with a quadratic dependence 
47r77/s = 3.64(r/ro - 1) + [T/Tof (blue dashed line) re- 
sembling the IQCD in quenched approximation in Ref. 
|46| as given also in [43|. We also consider a common 
case of a constant 77/s at its conjectured minimum value 
l/47r. 

In Fig. [21 we plot the results for the V2{pt) for the 
three different beam energies at RHIC and LHC at the 



same centrality 10 — 20% and for two different r]/s(T) as 
described in the following. On the left panel, the results 
are shown for 77/s = l/47r all over the evolution of the 
system. We see clearly that such a case would not predict 
an invariance of W2(pt) up to the LHC energy, but would 
generate a breaking up of about 20%. It is interesting 
also to notice that remaining in the regime of RHIC one 
would have indeed an approximate scaling of W2(pt) that 
however results to be misleading on a wider energy scale. 
This makes us to understand the importance of a wide 
beam energy scan (BES) up to LHC energy. 

In Fig. [2] (left panel) , we show by the dotted line the 
effect of modifying the initial pt distribution at LHC 
discarding the minijets. We clearly see that just this 
change can significantly affect the V2{pt) at least at pt > 
1.5 GeV. Therefore, behind the observed scaling there is 
an implicit role of the initial pT distribution that has to 
be correctly implemented including the non-equilibrium 
at increasing pT- Generally, the effect is that a stiffer 
distribution (like the one of mini-jets) produces a smaller 
V2{pt) even if the lifetime of the fireball and the assumed 
7;/s are unchanged. 

In Fig. [5] (right panel), we show the pattern of V2{pt) 
when an increase of rj/ s[T) is assumed in the cross-over 
region, solid or dashed line in Fig. [1] We label this case 
as Anrj/s = 1 + f.o. to emphasize that accounting for the 
increase of 77/s in the cross-over region naturally realize a 
freeze-out because it implies a smooth switching-off of the 
scattering cross section. For such a case, we see that the 
V2{pt) at different energies becomes more similar, even 
if still far from a scaling as observed experimentally. The 
main reason behind a more similar V2{pt) is that the fire- 
ball created at 62.4 GeV is more affected by the increase 
of 77/s in the cross-over region, while the system created 
at LHC is practically not affected at all. This of course 
makes the 772 (pt) more similar between RHIC and LHC, 
even if at 200 GeV it still remains larger because it is less 
affected by the increase of 'q/s{T) at low T compared to 
the 62.4 GeV case. 

The different impact oir]/s{T) on V2{pt) is determined 
by the different initial temperature and consequent life- 
time of the stage at T > Tc. In fact at RHIC energies 
such a lifetime is about 4 — 6 fm/c while at LHC about 
lOfm/c, in agreement with HBT results [4g]. Therefore 
at RHIC the elliptic flow has not enough time to fully 
develop in the QGP phase, while at LHC the lifetime 
is long enough to let the V2 develop almost completely 
in the QGP phase and the increase of rj/ s(T) at low T 
becomes irrelevant, as also found in Ref. 49]. 

From this reasoning one has the hint that an invariant 
V2{pt) can be caused by the specific T-dependence of 77/s 
that balances the suppression due to the viscosity above 
and below Tc where a minimum in 77/s should occur. We 
have considered the T-dependence of 77/s also at T > To 
in the QGP stage. In one case, we consider an rj/ s{T) 
rapidly (quadratically) increasing in the QGP phase, see 
the dashed blue line in Fig. [1] and in the other case a 
linear increase with T in agreement with IQCD data of 



Ref. [25|, corresponding to the red solid line in Fig. [T] 

As we can see comparing the right panel of Figs. [2] and 
the left one of Fig. [3]the rapidly increasing -q/ s{T) affects 
more the system created at LHC and this generates again 
a larger splitting of the V2{pt) among the different ener- 
gies. This again means that also a strong T-dependence 
in the QGP phase is in contrast with the observed V2{pt) 
invariance, essentially it would make the system at LHC 
too much viscous. Finally we consider Anrj/s = T/Tq, 
according to the solid red line in Fig. [TJ The results 
in the right panel of Fig. |3] are also compared with the 
experimental results for the V2 [4 1 ni easured at RHIC and 
LHC energy, data taken by |26l.l28||. We can see that in 
such case there is an almost perfect invariance of W2(pt) 
(within a 5%) in agreement with what is observed in 
the experimental data shown by the different symbols in 
Fig. [3] (right panel). The main effect is that with a mild 
increase of r]/s{T) at T > Tq the elliptic flow at LHC 
energies goes up reaching the higher V2{pt) obtained at 
lower energies. 

From a comparison with the first case considered in 
Fig. [5] (left panel), we understand that to have an in- 
variant V2{pt) it is essential also the rise of the r]/s{T) 
in the hadronic or cross-over region that significantly re- 
duces the V2 at low RHIC energy. Therefore only a fall 
and rise of the rj/ s{T) can account for a V2{pt) almost 
invariant going 62.4 GeV to 2.76 TeV. Furthermore, a 
comparison of the results for all the rj/ s{T) considered 
for RHIC at 200 AGeV shows that this is the case less 
affected by the T dependence of 77/s. Paradoxically this 
is the case more thoroughly studied by means of both hy- 
drodynamical and transport approaches till now. How- 
ever, we also notice that the impact of the T dependence 
of rj/ s on the V2(pt) is anyway quite weak and we have 
been able to educe the necessity of non vanishing T de- 
pendence only thanks to the experimental observation of 
the V2 scaling and exploiting a direct comparison in a 
quite wide range of colliding energy. Still one should no- 
tice that we have been able to discriminate a constant 
77/ s or a strong dependence like the quadratic one that 
for the maximum initial temperature at LHC would mean 
about a factor ten larger 77/s respect to the conjectured 
l/47r lower bound. The main reason is probably that 
even if the r]/ s is large at larger temperatures this is re- 
ally relevant only for the inner side of the fireball at the 
beginning of the expansion while most the elliptic flow is 
anyway formed later and more in the peripheral region 
of the fireball where the temperatures or energy densities 
are quite similar both as a function centralities or beam 
energy. However, we have shown that still a comparative 
analysis at different beam energies is able to reveal an im- 
portant information on the r]/s(T), namely the necessity 
of a "U" shape of rj/ s{T) with a minimum slightly above 
Tq as expected when the matter undergoes a cross-over 
\VA . [l9| . Therefore this finding provides a nice evidence 
for the phase transition of matter created in uRHIC's. 

Our result shows also that a BES with relativistic 
heavy-ion collisions allows to infer key properties of the 
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FIG. 3: As in Fig. [2]but for two different ri/s{T): a quadratic 
dependence on the left panel and a linear one on the right 
panel. In the right panel data for V2 [4] measured by STAR and 
ALICE collaborations [2B.l2a| are shown by different symbols 
as indicated in the legend. 



created matter that are not otherwise accessible and 
hence BES is far from being a mere repetition of a sim- 
ilar experiment. We notice that in the present study we 
find that the increase of r]/ s{T) below the QGP phase 
is consistent with previous studies at intermediate ener- 
gies which further supports the reliability of the study 
of nuclear matter through HIC's. In this respect an ex- 
pert reader could wonder why we studied the elliptic flow 
from BES starting from 62.4 GeV while data are avail- 
able from 7.7 GeV. We mention that this is due to the 
fact that at lower energy the baryon chemical potential 
jiB is no longer negligible, as pointed out also in [l9| . 
This would probably imply that a non-vanishing vector 
potential is acting as discussed in Ref. [50| . 

Our study is a seminal analysis of the information that 
we can obtain from the huge experimental efforts of the 
last decade, and one would expect that on a similar foot- 
ing more stringent constraints can be obtained from an 
analysis of higher harmonics like W3 . As this implies the 
development of initial state fluctuations in the transport 
approach and this is out of the reach of the present study, 
even if a work in such a direction is already under devel- 
opment. 

We find that for the different beam energies considered 
the suppression of the elliptic flow due to the viscosity 
of the medium has different damping coming from the 
hadronic or QGP phase depending on the average energy 
of the system. In particular, we observe that at LHC the 
elliptic flow is much less damped by the hadronic phase 
allowing a better study of the QGP properties. Moreover 
we have found that going from RHIC to LHC energies it is 
possible to have a nearly invariant V2{pt) only if the rj/s 
has a fall and rise with a minimum around the transition 
Tc ^ 165 MeV a behavior expected when there is a phase 
transition or a cross-over. 
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